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gas with Dirac spectrum 

S. M. Kukhtaruk, V. A. Kochelap 
Department of Theoretical Physics, V.E. Lashkaryov Institute of 
Semiconductor Physics NASU, Pr. Nauki fl, Kiev 03028, Ukraine 

Solving the initial value problem for semiclassical equations that describe two-dimensional elec¬ 
trons with the Dirac spectrum we found that collective excitations of the electrons are composed 
by a few distinct components of the oscillations. There always exist sustained plasma oscillations 
with well known plasmon frequency ijpi{k). Additionally there are oscillations with the ’carrier 
frequency’ uj = vpk slowly decaying in time according to a power law {vp and k are the Fermi 
velocity and wavevector). The reason for onset of these oscillations has a fundamental character 
related to branching of the polarization function of the Dirac electrons. A strongly anisotropic initial 
disturbance of the electron distribution generates additional component of undamped oscillations in 
the form of an electron unidirectional beam which are van Kampen’s modes in the Dirac plasma. 
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Polarization properties of electrons determine their re¬ 
sponse to an external electrical signal, define collective 
plasmon and plasmon-phonon modes, as well as many 
other phenomena involving interaction of the charges. 
A fundamental and applied aspects of collective exci¬ 
tations in graphene-like systems were discussed in de¬ 
tail in several reviewsic— . The polarizability of two- 
dimensional (2D) electrons with the linear Dirac spec¬ 
trum, e(p) = vfP, differs considerably from that of the 
electrons with a parabolic spectrum (e(p)ocp^) char¬ 
acteristic for conventional low-dimensional heterostruc¬ 
tures. Particularly, the absence of the spectrum cur¬ 
vature leads to less effective screening of the Coulomb 
potential^, gives rise to the divergence of the irre¬ 
ducible polarizability^— n(w, fc) oc Yj\JoP- — v'pk"^ at 
w -T ±VFk with u! and k are frequency and the 2D 
wavevector, respectively, vf is the Fermi velocity. Then, 
the square-root behavior of the polarizability means that 
n(a;, k) is a double-valued function on the complex oj- 
plane. For calculation of real characteristics of the phys¬ 
ical system, properties of the polarizability as an analyt¬ 
ical function on the complex w-plane is critically impor¬ 
tant. The semiclassical approach based on the analysis of 
the Boltzmann-Vlasov system of equations allows one to 
construct the principal branch of the polarizability func¬ 
tion n(a;, k) and obtain transparent and easy interpreted 
results. These results complement those calculated with 
the use of different quantum mechanical approaches in 
the long-wavelength limit^d^. 

Recently, new nanoscope techniques were proposed 
to launch the excitations by a sharp tip of atom force 
microscope and to monitor them by the scattering 
type near-held optical microscop o^^d^ . The techniques 
have enabled the experimental exploration of spatio- 
temporali^^— and time-resolved— dynamics of the col¬ 
lective excitations in the graphene and graphene-like sys¬ 
tems. Importantly, both the excitation amplitude and 
phase can be measured—. Similar studies have been 
performed by the use of the resonant antenna plas¬ 
mon launcher and near-held optical microscope^S. An¬ 


other promising method for plasmon investigation is 
time-resolved electrical measurement s^^ . Observed in 
these works macroscopic ehects - long-wavelength charge 
modes and local electric helds - can be described by using 
the semiclassical approach. 

In this Communication, we present the semiclassical 
analysis of longitudinal intraband excitations of 2D elec¬ 
tron gas with the Dirac spectrum. It is appropriate to 
note, that the transverse electric mode in graphene-like 
systems was analyzed in paper—. We assume the n- 
doped system and restrict ourselves to excitations with uj 
and k for which interband processes can be neglected at 
a given electron concentration, n, and ambient tempera¬ 
ture, T. Such approach is valid for hiv < cIjof), hk < pF 
at ksT <C e(j)F), with pF being the Fermi momentum 
and ks is the Boltzmann constant (detail discussion see, 
for example, in Ref. [l|). 

The Boltzmann-Vlasov system of equations consists of 
the collisionless transport equation 
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and the Poisson equation 
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[F{r,p,t) - Fo(p)] . 
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Here the Dirac spectrum is assumed for the electrons con¬ 
fined to the sheet at z = 0. The electron coordinate and 
momentum r, p are the 2D vectors in the x — y plane. 
F(r, p,t) is the electron distribution function, Fo(p) is 
that under equilibrium. $(r, z,t) is the self-consistent 
electrostatic potential. In Eq. e is the elementary 
charge, (5[z] is the Dirac delta-function, g is the degen¬ 
eracy factor of the electron band (for graphene g = A). 
The parameter k depends on dielectric environment. If 
the electron sheet is between two materials with the di¬ 
electric constants Ki, kh, in final formulae one shall use 
K= {ki+ Kh)l2,. 
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For the collisionless limit it is necessary ujt ^ 1, 
vpkr ^ 1, with r being a characteristic scattering time. 
For our purposes the equation m shall be linearized. 
We denote variation of the distribution function as 
J'(r,p,t) = F(r, p,t) - Fo(p) and set ^(r, p,t) = 
Jk(p, t) exp(ikr), and $(r, z,t) = $k(-z, t) ea;p(ikr). 

Now Eqs. o, and ([2]) take the form: 

+ WF—-Fk =-iek$k(0,t)^^ , (3) 
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Following the Landau approached, we consider the initial 
value problem by using the Laplace transform: 


OO 00+2(T 

/a;,k(p)=y-Fk(p,t)e*‘^*dt, Jk(p,t)= y/c.,k(p)e““*^, 

0 —oo-\-i(7 

(5) 

where cr > 0. Similarly, we define the transformation 
of the potential, (j)u,'k{z)- Now one can easily find the 
solution for /tj,k(p): 


/c.,k(p) 


.(5J'k(p) - *e(/ic^,k(0)kdFl)(p)/dp 
2_ 

w - •!;F(kp)/p 
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with (5J^k(p) being a given initial perturbation of the elec¬ 
tron distribution. The solution to Eq. 0, which decay¬ 
ing at z ± 00 , also can be easily found. The potential 
at z = 0 is: 


rh rn'i = [ dd'^P STi,{p) _ N{uj,'k) 

K fc A(u;, fc)y (27r?i)2 [a; — z;F(kp)/p] A(w,fc)’ 
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with A(w, fc) = 1 — 27re^ n(w, k)/nk and 
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It is clear that n(w, k) has the meaning of the polariz¬ 
ability obtained in the semiclassical limit. Fo(p) is the 
Eermi distribution, then 


n(a;, k) = 
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where Ep is the chemical potential. 

According to the definition of the Laplace transform 
©. in foregoing Eqs. 0-0 to = uj' + ioj" is the com¬ 
plex variable belonging to the upper half-plane, ui" > 0. 
Under the latter condition, the integral in Eq. 0 can be 
calculated as 


n(w, k) 


K K 
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The function dm is double-valued with two branch 
points, OJ = ivpk, as discussed in introduction. 

Consider an analytical continuation of n(a;, k) given 
by 0 to the lower half of the w-plane. We start with 
the use of (ttH) valid for all complex w-plane except the 
segment at the real axis —vpk <oj'< vpk. By choosing 
the cut along this segment, we select the principle branch 
of the function (HU. Defined by such a way, n(a;, k) 
is continuous when crossing the real axis at |w| > vpk, 
while it undergoes a jump in its value at the cut. 

As a result of this procedure, we obtain also the func¬ 
tion A(u;, k) as well-defined function in the complex oj- 
plane. It is easy to see that this function has simple zeros 
on the real axis 

OJ — '^ojpi , ojpi\k) — - -, (12) 

^(I + k/Kf - 1 

which are the exact solutions of the equation A(w, k) = 0. 
Eor what follows, it is important that ojpi(k) > vpk. 

Before to proceed with the farther analysis, let us es¬ 
timate the value K for the graphene. Assuming the 
graphene sheet over a AiC-substrate (k; pz 9.7, Kh = 1), 
setting the electron concentration n = 2 x I0^^cm“^ 
and temperature T = 300 K, we obtain Ep kj 0.16 eU, 
Ep/kpE Kj Q.12, pp/h = kp pz 2.5 x 10® cm“^, and 
K pz 3.9 X 10® cm~^. That is kp < K and for the semi¬ 
classical analysis we should use k < kp < K. For ex¬ 
ample, at k = 0.05 K we obtain ojpi pz 6.4 x 10^^ and 
Vpk PZ 2 X 10^® We refer to this set of parameters 
as K-environment-I. To estimate the criteria of validity 
of the collisionless approximation, one needs to know the 
characteristic scattering time, r. The ’intrinsic’ electron- 
electron scattering time is of the order of 10“^^ 3 ^. Op¬ 
tical phonon scattering can be neglected for the above 
accepted parameters. Acoustical phonon scattering time 
at T < 300 AT is estimated to be less than 2 x 10“^^ s^. 
Elastic scattering by imperfections is dominant. It can 
be determined via the phenomenological relationship for 
the transport time: r « pEp/ev'^ with p being the 
mobility. Assuming p = 2.5 x 10^ cm^/Us, we obtain 
T PZ A X 10“^® s and find that the semiclassical approxi¬ 
mation criteria are met: ojpiT Ri 26 ^ 1, vpkr « 8 1. 

Eor the graphene sheet in a high-K environment (for 
example, graphene in a solvent2£), the value K can be 
sufficiently less than kp. Then the case k> K may be 
actual. For example, at ki pz 9.7, Kh = 50, and k = K (k- 
environment-II) for the same electron concentration, we 
obtain K pz 7 x 10® cm“® and ojpi pz 8 x 10^® vpk pz 
7 X 10®® s“® (note, the accepted k is greater than for the 
K-environment-I). These estimates show that the excita¬ 
tions considered are rather of the THz diapason. 

Now, one can perform the inverse Laplace transform 
to find desirable functions in the time-domain. Below, 
we will concentrate on the potential at z = 0: 

oo+ 2 cr 

4>k(t) = $k(0,t)= y^,,k(0)e-“*^. (13) 
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According to Eq. 0, one can present (^aj,k(0) as a 
fraction, where the nominator, Af(a;,k), is determined 
by the initial perturbation, First, we will use 

the typical assumption^! that the initial perturbation is 
such that N{uj, k) has no poles in the w-plane. Then, the 
analytical properties of both A(w, k) and whole integrand 
in m allow one to deform the integration contour and 
calculate explicitly contributions of the residues related 
to the zeros of A(a;,k) and branch points: 


= $k (i) + ^k(i) 


= -*■ 


2 _ 2 p 

UJpi VpK 


3/2 


v'jrkK 
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k) exp [-iujpi t] - A(-Wpi, k) exp [iujpi t] 

duj 
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In the integral (fT6l) . contour C encloses the cut, as shown 
in the inset of Fig. 1(a). 

In Eq. (HH), the first term oscillating with time de¬ 
scribes the undamped plasmons excited by the initial 
perturbation of the distribution function. The latter 
specifies only the magnitude of the plasmon excitations, 
while their frequencies do not depend on the initial per¬ 
turbation. They are determined by the zeros of A(a;, k) 
given by Eq. For k <C AT, one obtains ’the square 

root law’—!!, ujpi VF\/kK/2; in the opposite limit, 
Wpi ^ vpk. 

The last term in Eq. (11^1) arises due to the existence of 
branch points in A(a;, fc). For electron systems with regu¬ 
lar quadratic energy dispersion, there is no branching and 
such contribution does not exist. The time dependence 
of this term is defined by the initial perturbation. Let us 
consider a few examples of different initial perturbations. 

First, assume that the initial perturbation is isotropic, 
(5Jk(p) = dTk{p), i.e., for perturbed electron gas the 
average momentum and velocity are zero. Then, 


0 

(17) 

where the square root should be defined as principal 
branch on the w-plane with the above discussed cut. 
Now, the contribution of the residues to the potential 
equals 




^ {k) — Vpk'^ 




cos[ujpi{k) t]. (18) 


The contour integral of Eq. ([TC| can be easily estimated 
for two limiting cases: 


{ k 

— Ji{vpkt)/vpkt, k 
Jo^vpkt), k ^ k , 


(19) 
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Figure 1. (color online). Normalized potential for: (a), 

(c) - K-environment-I and (b), (d) - K-environment-II; (a), (b) 
- isotropic initial perturbation; (c), (d) - strongly anisotropic 
initial perturbation, solid line - oo = tt/S, dotted line - oo = 0. 
Insets on (a), (c): cuts and integral paths on the u;-plane; inset 
of (b): for weak anisotropic perturbation; inset of (d): 

example of temporal pattern for ujpi « 2 uJvk- 


where Jo (a;), Ji(a:) are the Bessel functions. M t ^ 
l/{vFk)^ both expressions can be further simplified 




k cos[vFkt — 37r/4] 

2 I K {vFkt)^/'^ 

cos[vFkt — TT/4] 

[uFkiy/"^ ’ 


k<^K, 


k:^ K. 
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Thus, the initial isotropic perturbation of the electron 
distribution generates two different components of elec¬ 
trostatic potential oscillating in time and space. The first 
component is, obviously, sustained regular plasmon oscil¬ 
lations, excited by the initial perturbation. The second 
component corresponds to oscillations with the ’carrier’ 
frequency uj = VFk. The oscillations decay in time ac¬ 
cording to a power law. Mathematically, they arise due 
to the existence of the branch points of polarization func¬ 
tion ([5]). Recovering the space-dependent factor exp[ikr], 
one can see that 4)^ and components correspond to 
standing waves. The main panels of Figs. 1(a) and (b) 
show the normalized potential (p^{t) = 4)^/$° of these 
oscillations for two K-environments discussed above. 

Interestingly, the ratio of the magnitudes of the two 
components depends essentially on the wavevector. In¬ 
deed, for k <C K, and t —>• 0 we obtain the estimate 
ip9 Ri k/2K <C 1. From Fig. 1(b) we see that at k = K 
« 0.33. In the opposite case, fc AT, we obtain 
ip'^ Ri 1. That is, the regular plasmon component is 
preferable excited by long-wavelength perturbations. 

Then, consider an initial perturbation of the 
anisotropic form (5Jk(p) = '^•J'fe(p)pn/p, with n being a 
unit vector of a preferential direction. We denote the 
angle between n and k as ao- For this case, at < = 0 
the perturbed electron gas receives additional momen¬ 
tum and nonzero velocity, while the electron density is 
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unperturbed. The calculation gives 


N{u!, k) 




I cos ao 
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( 21 ) 


where $2 formally is given by the second relationship 
of Eq. (ITTI) . Calculating the time-dependent potential 
$k(i) with this fV(a;,k) function, one can use the ana¬ 
lytical continuation discussed above, as well as Eqs. (d 
- d, and the integration contour shown in the inset of 
Fig. 1(a). The residue contribution to the potential is 




= —2i cos ao 


v^kK^ 


sin[wpz(fc) t]. 
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For the contour integral of Eq. d, we present the results 
obtained for two limiting cases: 
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UJ 
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Ji{vpkt), 
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At t = 0, both contributions 4’^ and <1’^ equal zero. At 
t ^ l/vpk, we obtain: 





/ fcVcos [vpkt — ITm 

\K J {vpkty/'^ 
cos [upki — 37r/4] 


{vFkty/^ 


k 

k->K. 


(24) 

Thus, weak anisotropic initial perturbation also produces 
two component oscillations with purely imaginary mag¬ 
nitude, <1>^ and with frequencies w = ujpi{k) and 
oj = Vpk, respectively. The latter are decaying in time. 
At small t « 1 /vpk, these oscillations are quite differ¬ 
ent from those discussed above, since such perturbation 
does not generate directly a space charge and an electro¬ 
static potential. Both are developing with time. The nor¬ 
malized potential of these oscillations is 

illustrated in the inset of Fig. 1(b). The oscillations mag¬ 
nitude depends on the angle ao. The magnitude reach 
maxima for ao = 0, tt. If a = 7r/2, the perturbation does 
not excite the oscillations. 

Finally, consider the extremely anisotropic initial per¬ 
turbation: i5J^k(p) = ‘^pSTkip) S[a — ao], where a is the 
angle between p and k, ao determines the anisotropy 
direction. For the function A"(a;,k) we obtain 


N{uj, k) 


Lu — Vpk cos ao ’ 


(25) 


with still defined by Eq. (flT)) . The function (|25|) has 
a pole at the w'-axis: ui = uJvk{^) = vpk cosaQ . That 
implies that the accepted above way of the analytical con¬ 
tinuation of the function n(a;, k) is no longer applicable. 
Instead, one can use two cuts along semi-infinite lines in 
the lower part of the w-plane: uj = ±vpk + iut , uj < 0, 
as shown in the inset of Fig. 1(c). Now n(w, fc) is con¬ 
tinuous across the real axis everywhere excluding points 


LO = ±vpk and undergoes jumps at the semi-infinite cuts. 
The time-dependent total potential (IT^ consists of four 
contributions: 


<l>k(I) = <h«(t) + ch^Hi) + + $K^(t). (26) 


In contrast to the previous cases, every contribution has 
nonzero real and imaginary parts. The term $^(t) is 
the plasmon contribution determined by Eq. m with 
N{Lo,k) given by (|25|) . 


2 

uoh - v\k'^ 

3/2 

VpkK 


cos^ ao) 


x{ujpi cos[a;p;f] — ivpk cosao sin [ujpit]}. 


(27) 


The terms 4>^‘(t) and <i)^''(t), are the contour integrals 
around the left and right cuts, respectively (see Fig. 1(c)). 
Main panels of Figs. 1(c) and (d) display real parts of 
the total contribution -|- from the contour 

integrals. Latter can be found explicitly in the limiting 
case t :$> l/{vp k sin^ ao): 


4’k(0 


‘J’? 


sin^ao 


" k r 

— < cos [vpkt — 3 p/A] — i cos aoX 
K ( 

xsm[vpkt — Zp/4^l {vpkt)^^‘^,k<^K, 
^|cos [vpkt — p/A\ — i cosaox 
xsin [vpkt — p/4:]^^{vpkty^^,k:^K. 


The last term in 
at w = ujyKik)'- 


(28) 

is generated by the additional pole 




exp [-ioj^K{k)t] 

1 + K{i cosao/1 sinaoj -I- l)/fc’ 


(29) 


Thus, the strongly anisotropic initial perturbation pro¬ 
duces four components of the oscillations. As in previous 
cases, one of these components, ^^(f), is the sustained 
plasmon oscillations. Two other components, {t) and 
$^'■(1), are oscillations with ’carrier frequency’ w = vpk, 
they have different amplitudes and time dependences at 
finite t and both are similarly decaying in time. Restor¬ 
ing the space-dependent factor exp[ik.r\, one can see that 
4’^ and components correspond to traveling waves. 

Additional fourth component, corresponds to 

undamped oscillations with the frequency a;„ ^(k). This 
wave traveling along the wavevector k (or in opposite 
direction) can be interpreted as a particle beam mod¬ 
ulated in time and space and moving with the veloc¬ 
ity wf cosao. In physics of three-dimensional plasma 
with the regular dispersion of the electrons, such type of 
waves, with real w and fc, are known as the van Kampen 
modes^ii^. Particularly, van Kampen has showit^ that 
sharply anisotropic initial disturbance of the plasma gen¬ 
erates waves among which there is such a mono-energetic 
modulated wave. To distinguish the solution (1^^ from 
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other components of oscillations, we will designate it as 
a van Kampen mode. 

At a given fc, the van Kampen mode depends on two 
parameters, and ao, which characterize the initial per¬ 
turbation. When the angle of the anisotropy direction, 
ao, changes from 7r/2 to 0, the frequency of this mode 
increases from 0 io vpk. At that the magnitude of the 
wave varies from infinity (ao = 7r/2) to zero (ao 0). 
In the latter case, when the pole vpk, the 

contribution sharply increases as illustrated by numerical 
calculations presented in Figs. 1(c), (d). Interestingly, in 
the Dirac plasma, the modulated electron beam corre¬ 
sponding to the van Kampen mode is composed by uni¬ 
directional but not necessary mono-energetic electrons, 
as it is seen from the accepted above form of the initial 
perturbation of the electron distribution. 

Note, the region on the ui-k plane, where the van Kam¬ 
pen modes can be excited, coincides formally with the 
region, for which the quantum theory predicts the exis¬ 
tence of the electron-hole pairs at low temperatureii^^— . 
As contrasted with such one-particle excitations, the van 
Kampen modes are collective charge excitations. 

At a given k, all obtained oscillating components are 
spatially coherent. Their superposition (1261) may demon¬ 
strate a complex temporal behavior, for example, a beat¬ 
ing effect between the plasmon waves and the van Kam¬ 
pen modes, as shown in the inset of Fig. 1(d). To char¬ 
acterize a complex temporal signal, it is useful to apply 
the Fourier analysis. We performed such an analysis for 
obtained oscillations (1261) and found the amplitude and 
phase of this signal as functions of the frequency. The 
amplitude shows very narrow peaks at frequencies of the 
regular plasmons and the van Kampen modes, while re¬ 
spective phases behave with sharp jumps, as it is typical 


near resonance. At the frequency uj = vpk, correspond¬ 
ing to the component of oscillations, the amplitude 
has a smeared spike, however the phase appears clear 
characteristic kink. These features of the amplitude and 
phase of the collective excitations can be verified by con¬ 
temporary measurement techniques. 

In conclusion, we have used the semiclassical approach 
to obtain transparent results on the dynamics of collec¬ 
tive excitations in Dirac 2D electron gas. By solving 
the initial value problem for the system of Boltzmann- 
Vlasov and Poisson equations with different forms of ini¬ 
tial disturbances of the distribution function and charge, 
we found that collective excitations of the electrons with 
the Dirac spectrum are composed by a few distinct com¬ 
ponents of the oscillations. Among these, there are al¬ 
ways well known sustained plasmon oscillations with the 
frequency given by Eq. dUD. Also, there exist new type 
of oscillations with the carrier frequency uj = vpk' they 
are of a transient character, decaying in time according 
to a power law. Mathematically, the latter component of 
the oscillations arises due to the branching feature of the 
polarizability function. Finally, a strongly anisotropic 
initial disturbance generates another component of un¬ 
damped oscillations of frequency uj = VFkcosao, with 
ao being the angle between the wavevector k and the 
anisotropy direction. We interpreted these undamped os¬ 
cillations in the form of an electron unidirectional beam 
and attendant electrostatic potential both modulated in 
time and space, as van Kampen’s mode^l in the plasma 
of the Dirac electrons. 

The results obtained demonstrate that the semiclassi¬ 
cal approach is adequate to describe the dynamics of col¬ 
lective excitations of THz spectral diapason in the Dirac 
electron gas. 
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